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In the recently discovered Co-oxide superconductor Na^CoC^ • {/H2O, the edge-shared C0O6 
octahedra are trigonally contracted along the c-axis in the Co02-plane. We study how this CoOe 
distortion affects the magnetic properties and superconductivity in this compound by analyzing 
the multiorbital Hubbard model using the fluctuation-exchange approximation. It is shown that 
through generating the trigonal crystal field, the distortion pushes the Co e' s bands up and 
consequently gives rise to the hole-pocket Fermi surfaces, which have been predicted in the band 
calculations. As the distortion increases, the hole pockets are enlarged and the ferromagnetic 
fluctuation as well as the pairing instability increases, which is in good agreement with recent 
NQR results. 

KEYWORDS: Co-oxide superconductor, hole-pocket Fermi surfaces, ferromagnetic fluctuation, triplet pairing 
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The recently discovered Co-oxide superconductor 
Naa;Co02 • 2/H2O 1 ) has attracted great interest. One 
of the interesting points is that the unconventional su- 
perconductivity due to an electron-correlation mecha- 
nism, in particular, the spin-fluctuation-mediated mech- 
anism, is expected in this material. Although there are 
some controversies, NMR/NQR 2-4 ' and /i,SR 5,6 ' mea- 
surements have shown evidence for non-s-wave pairing. 
Various experiments on the nonhydrate compounds show 
characteristic behaviors of strongly correlated electron 
systems. 7 ' 8 ^ Thus, the magnetic properties have been in- 
tensively studied since we recognized that understanding 
them is essential for elucidating the pairing mechanism. 
A ferromagnetic (FM) spin fluctuation was predicted 
in the LSDA calculation 9 ' and was claimeded based on 
NQR experiments by several groups. 2 ' 3 ' 10 ' In contrast, 
some groups observed Curie- Weiss behavior of 1/Ti 4 ' 
and a temperature-independent Knight shift 11 ', which 
is evidence against the FM fluctuation. The magnetic 
properties of this compound are still controversial. 

Quite recently, NQR measurements performed by 
Ihara et al. showed that there is a correlation between 
the magnetic fluctuation and NQR frequency vq aris- 
ing from the ±5/2 <-> ±7/2 transition. 12 ' 13 ' They found 
that the higher-z^Q sample has a stronger magnetic fluc- 
tuation and a higher superconducting transition temper- 
ature T c . In addition, a weak magnetic order was found 
in a sample with the highest vq. They proposed that 
the magnetic fluctuation and T c have a correlation with 
the distortion of the CoOg octahedra along the c-axis 
since vq is expected to scale with the CoOg distortion 
from the cubic symmetry. They also proposed that su- 
perconductivity occurs in the vicinity of magnetic order, 
suggesting a spin-fluctuation-mediated pairing. A rela- 
tionship between T c and the c-axis parameter has also 
been pointed out by several groups. 14 ' 15 ' 

Another interesting point is orbital degeneracy. This 
compound has quasi-two-dimensional C0O2 planes with 



edge-shared networks of CoOg octahedra, in which the 
Co ions form a triangular lattice. Band calculations 
showed that the Fermi surface is constructed from two 
bands consisting of three Co t 2g orbitals. 16 ' Concretely, 
it was shown that among these three t 2g orbitals, the 
singlet aig orbital forms a large cylindrical Fermi sur- 
face around the T-point and six small hole pockets are 
constructed from the doublet e' g orbitals near the K- 
points. Considering these aspects, we have recently pro- 
posed that the Naa;Co02 system should be described by 
a multiorbital model 17-19 '; most of the previous works 
are based on single-band models. 20-24 ' In the studies 
reported in refs. 17 and 18, we analyzed a multior- 
bital Hubbard model by using the fluctuation-exchange 
(FLEX) approximation and perturbation theories, and 
found several important aspects which had not been 
expected in the single-band theories. These analyses 
showed that when the system has the six hole pock- 
ets predicted in the band calculations, the interorbital 
Hund's-rule coupling induces FM fluctuation and this 
FM fluctuation favors the spin-triplet pairings with 
fy( y 2_ 3x 2)-wa,ve and p-wave symmetries. We discussed 
that for the FM fluctuation as well as the triplet super- 
conductivity, the presence of e' g hole pockets is crucially 
important, which was also pointed out by Kuroki et al. 23 ' 
and by Yata et al. 25 '. 

In this letter, we study the results of the NQR mea- 
surements preformed by Ihara et al. 12,13 ' by extending 
our previous study on the multiorbital Hubbard model 
in the FLEX approximation. We show that the results 
of Ihara et al. can be clearly explained if we consider 
the hole-pocket Fermi surfaces as well as the Co t 2g or- 
bital degrees of freedom. In the C0O2 plane, the C0O6 
octahedron is trigonally contracted along the c-axis as 
shown in Fig. 1(a). As a result, the O ions generate a 
trigonal crystal field (TCF) on the Co site, which lifts 
the local Co t 2g degeneracy into the lower singlet a\ g 
and higher doublet e' g orbitals (see Fig. 1(b)). First, we 



2 



Masahito Mochizuki, Youichi Yanase and Masao OGATA 



(a) 



(b) 




C O ion 



-eg 
A 

■ a\% 

cubic trigonal 



(°) 2 
o 

-2 
-4 



hole pocket 









A=0.8 







K M 



Fig. 1. (a) Trigonal distortion of C0O6 octahedron and a\ 9 or- 
bital. Arrows indicate shifts of O ions, (b) a\ 3 -e' g splitting of t2 g 
degeneracy due to trigonal crystal field, (c) Band dispersions cal- 
culated from noninteracting Hamiltonian H^ n + H CTy . . The pa- 
rameters are chosen to be ti = 0.38, = 0.20, t$ = 1, £4 = 0.12, 
t 6 = -0.15, te = -0.06, t 7 = 0.12, t 8 = 0.12, = -0.30, 
and A = 0.8. The electron number is 5.55 per Co site. The 
linewidth of each band in (a) is proportional to the amount of 
the ai 9 -orbital component. 



show that when the a\ g -e! g splitting becomes larger with 
increasing CoC>6 distortion, the hole pockets appear and 
become larger near the K-points. As the hole pockets be- 
come larger, the FM spin fluctuation becomes stronger 
and the pairing instability is enhanced, while they are 
markedly weak without the hole pockets. The successful 
reproduction of the experimental results seems to sup- 
port the existence of hole-pocket Fermi surfaces and the 
pairing mechanism proposed in the previous papers. 17 ' 18 - ) 
The multiorbital Hubbard Hamiltonian is given 



with H, 



cry. 



by H — H CTy . + iJkin. + -ffint. 

En J tt v-> j-mn J\ J 

i,m,n,a ±J mn u i m< y u 'in<Tt -"kin. — Z-,i,j,m,n,<7 °ij ""ima^jna 

M 2i,m,<T A-ma^ima , and H int . = H\J + HjJ' + H j n + Hj' , 

where di ma (d\ ma ) is the annihilation (creation) operator 
of an electron with spin c(=T, [) in the orbital m on the 
Co site i. Here, m runs over the xy, yz, and zx orbitals. 

The first term H cry , expresses the TCF from the O 
ions. The matrix element D mn is given by ^{5 mn — 1), 
which gives the a\ g -e' g splitting of A. The second term 
i?kin. is a tight-binding Hamiltonian, where the hopping 
integral i expresses the electron hopping between the 
m orbital on the i site and the n orbital on the j site. 
We choose the lattice constant as a unit length and de- 
note the unit vector as d{=(-\/3/2,— 1/2) and O2=(0,l). 
We consider up to the third nearest-neighbor hoppings 
between orbitals. In the fc-space, -f/kin. is rewritten as 



-ffkin. — 



with 



,77 _ 



k,m,7i,er c k a kmcr 



A 4 Ek.m.a Ai-ma^ma 



2*i cos kji + 2t 2 [cos fc 77 + cos(fc^ + fc 77 )] + 
2U [cos(2/c2 7 + k^ 1 ) + cos(/c2; 7 - ^ 7 )]+2i 5 cos2fc2 7 , and 
e JT' = 2i 3 cos/c 77 ' + 2i 6 cos2fc 77 ' + 2i 7 cos(fc2 7 '+2/c 77 ') + 
2i 8 cos(fc2; 7 '-fc 77 ') + 2i 9 c os(2A:2 7 ' + ^ 7 '). Here, 7 and 7' 
represent xy, yz and zx orbitals and k* v ' xy = k* v,zx = 

ti, k™ = k x b v > zx = k 2 , ky*> yz = k x a y^- = fc 2 , k%*' vz = 

(fci + *a), k z a x ' zx = fcf ' zx = -{ki + k 2 ), and 
fci, respectively. k\ = \/3/2k x — l/2k y 
and k 2 = k y are the components of the wave vector k of 
the triangular lattice spanned by d{ and 02, respectively. 
In Fig. 1(c), we show the band dispersions obtained from 
#kin. +H cxy . with ti = 0.38, t 2 = 0.20, t 3 = l,t 4 = 0.12, 
is = -0.15, i 6 = -0.06, i 7 = 0.12, i 8 = 0.12, i 9 = -0.30, 
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and A = 0.8. In the following, we use these parameters 
and *3 = 1 as the energy unit. The LDA results are well 
reproduced, particularly near the Fermi level. 17 - 16 ) The 
linewidth of each band is proportional to the amount 
of the ai g -orbital component, which shows that the hole 
pockets have an e' ff -orbital character. Here, we note that 
the dispersions far below the Fermi level are relatively 
different from those of LDA. This may imply that it is 
necessary to consider the electron hoppings beyond the 
third nearest neighbors or the O 2p orbital components 
to reproduce the precise structure of the LDA result. 
However, we have confirmed that the results in this pa- 
per do not depend on details of the band structure far 
below the Fermi level. 

The last term in H, ffint., represents the on-site d-d 
Coulomb interactions, where Hjj and Hjji are the in- 
traorbital and interorbital Coulomb interactions, respec- 
tively, and Hj m and Hj> are the Hund's-rule coupling 
and the pair-hopping interactions, respectively. These 
interactions are expressed using Kanamori parameters, 
U, U' , Jh, and J', which satisfy the relations U' — 
U — 2Jh and Jh = J'. The value of U has been es- 
timated as 3-5.5 eV by photoemission spectroscopy, 28 -* 
and the value of Jh for the Co 3+ ion is 0.84 eV. Thus, 
the ratio Jn/U, which gives the strength of Hund's-rule 
coupling, is 0.15-0.28 in this compound. 

We analyze this multiorbital Hubbard model by ex- 
tending the FLEX approximation to the multiorbital 
case. 17 ' 29 -* In the present three-orbital case, the Green's 
function G, the noninteracting Green's function and 
the self-energy E are expressed in 3 x 3-matrix form cor- 
responding to the xy, yz, and zx orbitals. 17 ' Calculations 
are numerically carried out with 64x64 fc-meshes in the 
first Brillouin zone and 1024 Matsubara frequencies. The 
value of U is fixed at 8.0, and temperature T is fixed at 
0.02. 

First, we discuss the change of the band structure when 
the trigonal CoOe distortion is increased. Since the TCF 
pushes the e' g level up, we naively expect that the hole- 
pocket Fermi surfaces would appear as TCF splitting A 
is increased. In Fig. 2, we show the band dispersions (a) 
and Fermi surfaces (b) in the case of Ja/U =0.18 and 
T = 0.02 for various A values obtained in the FLEX 
approximation. They are calculated by diagonalizing a 
matrix whose component is e£? n + E m „(k, ui = 0). Here, 
the self-energy is given by 



q \1V 



M)G^{k-q) (1) 



with effective interaction V^ m ,vn{(l) in RPA-type bubble 
and ladder diagrams. 17 ^ Figure 2 actually shows that the 
e' g bands are gradually pushed up with increasing A. As 
a result, when A >0.8, the e' g bands intersect the Fermi 
level and form the six hole pockets near the K-points. 
With further increasing A, the areas of the hole pockets 
become larger. 

We next discuss the correlation between the spin fluctu- 
ation and the magnitude of the trigonal distortion. Fig- 
ure 3 shows the spin susceptibility \ s (q) = X s (q> w = 0) 
for various values of A at Jh/£7=0.18 and T =0.02. 
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Fig. 2. (a) Band dispersions and (b) Fermi surfaces for various 
A. The parameters are fixed at U = 8.0, Jn/U = 0.18 and 
T = 0.02. 
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Fig. 3. (a) Spin susceptibility x s (q) f° r 0< A <1.2 and (b) 
X s (q) for 0.6< A <1.8. The parameters are fixed at U = 8.0, 
J H /U=0.18, and T =0.02. 



Here, x'C?) is given by x s (?) = Em„Xmm,n«(?) with 
X s (?) = [/ — X°(9)^ S ] _1 X°(9) an d C7 S is a 6x6 matrix 
representing the interactions U, U', Jjj, and J'. 17 ) Fig- 
ure 3(a) shows that the nature of spin fluctuation is 
completely different between the case with hole pockets 
(A > 0.8) and the case without hole pockets (A < 0.6). 
X s (q) for A > 0.8 has an enhanced peak structure at the 
T-point indicating the dominant FM fluctuation. This 
FM fluctuation rapidly increases with increasing A as 
shown in Fig. 3(b). On the other hand, without hole 
pockets, x s (q) does not have any remarkable structure 
except for very tiny peaks caused by the weak nesting 
of the dig Fermi surface. Since the Hund's-rule coupling 



Fig. 4. Eigenvalues of Eliashberg 
and p-wave pairing at U = 8.0, 
functions of TCF splitting A. 



equation for /j / (j / 2_ 3a .2)-wave 
J H /(/=0.20, and T =0.02 as 



only works between two electrons on different orbitals, 
the coupling cannot generate the FM fluctuation when 
the system only has the a\ g Fermi surface. These results 
suggest that the existence of hole pockets is essential for 
the appearance of the FM fluctuation. 

To discuss the nature of superconductivity, we solve 
the following Eliashberg equation: 

T 
N 



\4>r, 



= EEC 

a/3 fif 



^(Q)4>i^(k - q) 



x Gaji (k - q) Gj3v{q - k). 



(2) 



The singlet and triplet pairing interactions T s and 
T* have 9x9- matrix forms as T s (q) — §f7 s x s (g)?7 s — 
\U c x c {q)U c + \{U S + U c ) and f*^) = -±U s x s (q)U s - 
\U c x c {q)U c + \{U S + U c )- Here, X c {<l) = [I + 
X (g , )C/ c ]~ 1 x (g) and U c is a 6x6 matrix representing 
the on-site Coulomb interactions. 17 - ) The eigenvalue A is 
a measure of the dominant pairing symmetry, and be- 
comes unity at T = T c . 

In Figure 4, we plot the eigenvalue of the Eliashberg 
equation (A) for both f y ( y 2_ 3x 2yw&ve and p-wave states 
as functions of the TCF splitting A. As discussed previ- 
ously, 17, 18 ) the triplet f y ( y 2_ 3x 2ywa,ve or p-wave pairing, 
which has the large gap amplitude on the e' g hole pock- 
ets, has the largest A and the other pairing instabilities 
are markedly weak. We can see that A for the fy(y^-^x 2 )- 
wave state and that for the p-wave state are markedly 
small in the region A < 0.6, which indicates that with- 
out the hole pockets, the pairing instability is strongly 
suppressed. From A ~ 0.7, A rapidly increases as A 
increases. For A larger than ~1.4, the system exhibits 
FM ordering. Here, the FM phase transition is identi- 
fied by Min{det[/-x°(g)^ s ]} = 2 x 10" 3 . This behavior 
is quite consistent with the NQR result, which suggests 
the increase of T c with increasing CoOg distortion and is 
consistent with the finding of weak magnetic ordering in 
a sample with more significant distortion. 12 ' 13 ) 

In summary, motivated by the recent NQR experi- 
ments, 12,13 ) we studied how the C0O6 distortion along 
the c-axis affects the magnetic fluctuation and the pair- 
ing instability in Na^CoC^ • 2/H2O on the basis of the 
FLEX analysis of the multiorbital Hubbard model. We 
showed that through generating TCF, the CoC>6 distor- 
tion gives rise to the hole-pocket Fermi surfaces predicted 
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in the LDA calculation by pushing the e' g bands up. The 
enhanced FM spin fluctuation is realized when the sys- 
tem has the hole pockets, and the magnetic and pair- 
ing instabilities are enhanced as the distortion increases. 
These results are in good agreement with the experi- 
ments and, thus, seem to support the validity of our 
theoretical proposal about the triplet pairing being me- 
diated by the Hund's-rule-coupling-induced FM fluctu- 
ation, which assumed the existence of the hole-pocket 
Fermi surfaces. Further, the results suggest a possible 
controllability of T c and magnetism in this compound 
via the application of pressure along the c-axis. The 
pressure may possibly induce superconductivity even in 
the nonhydrate systems. 

Note that in angle-resolved photoemission spec- 
troscopy (ARPES) studies on nonhydrate compounds, 
the hole-pocket Fermi surfaces have not been observed 
thus far. 8 ' 26 - 27 ) This may due to the surface sensitivity 
of the CoC>6 distortion. Our calculations show that the 
trigonal crystal fields due to the contracted CoOg octa- 
hedra generate the hole pockets by pushing the e' g bands 
up. The contraction may be relaxed at the surface re- 
sulting in the vanishing of hole pockets. Although the 
existence pf hole pockets has not been confirmed exper- 
imentally at present, our results strongly support their 
presence. 

Within our model, the triplet f y ( y 2_ 3x 2ywa,ve and p- 
wave pairings are favorable. To elucidate the pairing 
state, several groups performed Knight-shift measure- 
ments by NMR 2,10 ' 30 ' 31 ) and fiSR 5 \ Some of these 
groups reported a constant Knight shift below T c . 2 ' 5 ' 31 ) 
This could be evidence for the triplet pairing. On 
the other hand, other groups reported a decreasing 
Knight shift below T c , which suggests the singlet pair- 
ing or triplet pairing whose <i-vector is fixed in the ab- 
plane. 10 ' 30 ) These two results seem to be contradictory 
at first sight. However, they may be understood from dif- 
ference in the experimental conditions. In fact, the for- 
mer measurements were performed under relatively high 
magnetic fields, while the latter ones were performed un- 
der rather low fields. These two results may suggest 
the possible existence of multiple superconducting phases 
under the magnetic field, and the d-vector may be di- 
rected along the c-axis in the high-field phase, while in 
the low-field phase it may be directed in the a6-plane. 
Quite recently, it was reported that the Knight shift de- 
creases even when the magnetic field is parallel to the c- 
axis. 32 ) . However, the decrease is much smaller than that 
expected from the c-axis magnetic susceptibility. Com- 
plete understanding of the NMR results requires further 
theoretical studies on the direction of the (i-vector by 
considering the effects of the spin-orbit interaction and 
magnetic fields. 
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